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Abstract 

Results of extensive Monte-Carlo simulations that investigate the out-of- 
equilibrium dynamics of the one-dimensional Ising spin glass model with a 
Gaussian bond-distribution are presented. At low enough temperatures a 
typical (interrupted) aging scenario is established as in two- and similar to 
three-dimensional spin glass models. Since the underlying mechanism is a 
slow domain-growth we study in detail spatial correlations and the time- 
dependence of the domain- as well as kink-statistics. We find that all cor- 
relation functions in time and in space as well as the domain-size probability 

distribution function obey simple scaling laws. 
75.10N, 75.50L, 75.40G. 
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I. INTRODUCTION 



Non-equilibrium dynamics in real spin glasses, also known as aging, has been a major 
focus of research interest of experimentalists since many years Jl|-§||. Quite recently, also 
theoreticians began to shift their attention from the equilibrium statistical mechanics 0] 
to the out-of-equilibrium dynamics of various spin glass models [|5|-|T2]j. It has been argued 
that aging is a characteristic feature of the low temperature dynamics of spin glasses and 
it has already been established in experiments |T3"|,H| and numerical simulations [§,T3] in- 
vestigating two-dimensional spin glasses (which do not posses a finite temperature phase 
transition) that a true spin glass phase is not a necessary prerequisite for it to be observable 
on suitable time-scales. The only differences remaining is then that aging is interrupted at 
some maximal waiting time, which is equivalent to the finite, but large, equilibration time 
at this temperature (see also @ and 0). 

In this paper we will demonstrate that even in a disordered but non-frustrated system 
a very similar scenario will occur. The model we consider here is the one- dimensional Ising 
spin glass with a Gaussian bond-distribution, which can be mapped, simply via a gauge- 
transformation of the spins, onto a random, ferromagnetic Ising chain. Thus the ground 
state is fully magnetized and the dynamics at low temperatures is characterized by (fer- 
romagnetic) domain-growth, which is drastically slowed down by the presence of a huge 



number of metastable states [|16f fL8|j. An important ingredience, however, is a continuous, 
unbounded bond distribution, although some of the features to be presented can be antic- 
ipated already in the one- dimensional ± J Ising spin glass model. The latter is isomorphic 
to the ferromagnetic Ising chain, whose non-equilibrium dynamics is exactly solvable and 
has been studied extensively [|l9H2"2|j. Only for a continuous bond distribution one has a 
broad distribution of (free) energy barriers, which are responsible for the broad spectrum of 
relaxation times. 

The paper is organized as follows: The next section reports known exact results for 
the autocorrelation function and the two-point correlation function of ±J Ising chain. In 
section 3 we present results of Monte-Carlo simulations of the random bond Ising chain, 
which cover correlations in time and in space. A scaling analysis along the lines sketched in 
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section 2 is performed and compared with the results of a recently proposed one-dimensional 
domain-growth model Section 4 is devoted to a domain-size and -kink statistics that 



we obtained from our numerical simulations and here also a scaling analysis is performed. 



Section 5 discusses the concept of an overlap length |20| , |2"5|| in the random bond Ising chain 
and summarizes our findings. 

II. THE ±J CHAIN 

In this section we recapitulate what is known about the non-equilibrium dynamics of the 
±J chain, which can be solved exactly |T9|-|23|. On one hand we introduce in this way the 



quantities that we intend to study for the random bond Ising chain and on the other we shall 
see what kind scaling behavior can be expected for these quantities. The one-dimensional 
± J Ising spin glass model (which is no spin glass at all) is defined by the Hamiltonian 

N 

H = —^JiSiSi+i , (1) 
i=i 

where Si — ±1 are Ising spins and the bonds Jj are chosen at random from a binary 
distribution 

P(Ji) = +6{Ji + l)], (2) 

however, this quenched disorder can be removed (up to a negligible boundary term) via the 
Gauge-transformation S' t = Si 111=1 sig n Jk- Thus ([I]) is equivalent to the simple ferromag- 
netic spin chain. All results reported below are for N — > oo. The dynamics is the usual 
Glauber-dynamics fj~9 |, where each spin is flipped with a probability 



w(Si -> -Si) = -[1 - tanh (3(Ji-iSi-i + JiS i+1 ) ] . (3) 

The corresponding Master equation leads to a time dependent probability distribution for 
spin configurations with which expectation values are calculated (the initial distribution is 
uniform, i.e. the at t = all spin are chosen at random). Such a quantity for instance is the 



autocorrelation function [fL9 l 



C{t,t w ) = (Si{t w )Si{t + t w )) , (4) 



which will be studied numerically for the random bond chain in the next section. At T = 
one obtains [ |T9| , [2T 



2 r in - 1/2 

C{t, t w ) = -arc sin \l + -— \ . (5) 



7T L 2t u 

Note that C(t,t w ) is a function of t/t w only and such a scaling behavior has been reported 
for the same quantity in various other spin glass models |^,[7|,p|, |II| , |l5| . For illustration and 
comparison with these models we have depicted in figure p] a plot of C(t,t w ) for various 
waiting times. Furthermore one gets for short and for long times 

{ 1 - | (t/t w y/ 2 for t « t w , 
C(t,t w ) « \ (6) 

[ ^ (t/t^ 1 / 2 for t » . 

For T > the time-dependence of C(t,t w ) is determined by the correlation length (which is 
infinite for T = 0), which we consider now. Spatial correlations are defined by 

G T (r,t w ) = (S l (t w )S l+r (t w )} . (7) 

For t w — > oo one has 

lim G T (r,t w ) = G c ^{r) = {tanh(/5J)} r (8) 

which defines the equilibrium correlation length £<r,eq via £^ eq = — log{tanh(/3J)} fl9f . It 



can be shown |2T] that for 1 <C t w -C 



.cq 



where the characteristic length scale grows like 

6r(*«0 oc (10) 



until it saturates for t w — > £f. and one recovers GT(r,t w ) — > G^(r). Similarly, for 1 <C 
"C ^T,eq one would expect instead of (|5|) 

Cr(M » )=ST (?pk)) ■ 

where the characteristic time scale grows like r(T,t w ) ~ £ M for small temperatures, but 
saturates at £| q (T). All the above results are very reminiscent of characteristic aging phe- 
nomena observed in other spin glass models. The main difference to these much more 



complex models lies in the waiting time dependence of the correlation length £r{tw), which 
predicts a rather fast domain growth, which makes it hard to observe them on macroscopic 
time scales at non-zero temperatures. Furthermore, a significantly different scaling behavior 
for e.g. C(t,t w ) will be present in systems that are characterized by activated dynamics and 



logarithmic domain growth |2~5*|. 



III. THE RANDOM BOND ISING CHAIN 

The random bond Ising chain is defined in the same as the ±J Ising chain discussed 
in the last section up to a different bond-distribution: Instead of (0) we consider here the 
Gaussian distribution of the random bonds 

P(J 4 ) = -L exp ( -£) , (12) 
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(see e.g. |23] for the corresponding spin glass model in higher dimensions). As before one 
can remove all minus-signs in the bonds by a gauge-transformation that leads to a purely 
ferromagnetic bond distribution P(Ji) = 2P(Ji)9(Ji). 

We performed Monte-Carlo simulations of the single-spin-flip Glauber dynamics @ with 
typical chain lengths of iV = 10000 (imposing periodic boundary conditions) and averaged 
over 128 to 512 samples, i.e. realizations of the disorder. Since we performed our simula- 
tions on a Parsytec GCell024 parallel computer, using up to 256 different processors. We 
constrained the system size to fit one system into the memory of one T805-transputer in 
order to avoid communication delays. Alternatively we also could have made a few runs for 
chain-lengths of 10 6 -10 7 and distribute chain- segments over several processors, however, by 
varying the system size, we checked that our date were not spoiled by finite-size effects. 

The correlation functions (^) and (|7|) have to be averaged over the quenched disorder. In 
addition, in order to obtain better statistics we calculated instead of single-site/single-time 
quantities, appropriate space and time averages as 

^ N i t+At 



(13) 

i N | t+At 

G ^ U = x E E [SiV + t w )s i+r (f + o] av , (i4) 



with suitable values for the time-window At (1 < At -C t) and [• • -] av indicating the disorder 
average. If taken literally the expectation values in (|J) and (|7|) should be obtained via 
different realizations of the thermal noise. However, we checked that in using (|IB,rj) we get 
the same results. 

In figure |2] we show our results for C(t,t w ) at different temperatures. The picture is 
equivalent to that obtained for the same quantity in the two-dimensional Ising spin-glass 
|15| and, at low temperatures, reminiscent of the three-dimensional spin glass model @ 



and the SK-model within the spin glass phase ]TT[. This is what we would like to call 



interrupted aging as been observed experimentally in charge-density wave systems ||. At 
higher temperatures the curves collapse for larger waiting times t w , which means that the 
system is equilibrated. The waiting time at which this happens is the equilibration time r eq . 
At lower temperatures this time scale is out of reach for the amount of CPU time that is 
available to us, hence all curves seem to be shifted logarithmically on the time axes by t w . 



Similar to what has been observed in other spin glass models ||E%[7||^ Jliyi5|] we established 
the scaling 

A scaling plot of this kind is depicted in figure where the insert shows the characteristic 
time-scale r(T, t w ). For small waiting times at low temperatures it is, as expected, r(T, t w ) « 
t w , and it saturates for increasing t w at r(T,t w ) = r eq (T). 

Next we turn our attention to the spatial two-point correlation function (|7|). The (aver- 
aged) equilibrium correlation function can easily be determined exactly as 

G cq (r) = [tanh(/3J)L (16) 

(for a discussion about the differences between averaged and typical or most probable values 



of the spatial correlations in this system see []26[]) which yields an equilibrium correlation 
length 

C q 1 (T) = -log[tanh(/3J)] av ocT for T « 1 . (17) 

The non-equilibrium spatial correlation function GT{r,t w ), defined in (|7|) respectively (|T4|). 
will approach (|T6|) at waiting times t w larger than the equilibration time r cq (T). This is 



shown in figure |4] for two different temperatures, where at T = 0.4 one observes indeed this 
data-collapse. The curves are approximate straight lines in a linear- log plot, which means 
that the spatial correlations decay exponentially with a characteristic length scale £r(^)- 
This length is roughly equal to the length scale <TT{t w ) that yields a good data-collapse for 
a scaling plot according to equation ([5J), i.e. 

G(r,t w ) =g(r/a(t w )) , (18) 

which is shown in figure [5| The smoothest curve for the waiting-time dependent length scale 
can be obtained via the introduction of the (effective) correlation length 

POO 

£t(U = / drG T (r,t w ) , (19) 
J o 

note, however, that for any t w and T we have £t(£«>) ~ ctt(^to) ~ iritw)- The correlation 
length £t(£«>) is shown in the insert of figure |5| within a log-log plot, where one observes that 
this length scale grows much slower than in the ±J chain, cf. (0)- Similar to what has been 



found in the two-dimensional spin glass model [TjJ it seems that (see the curve for T = 0.2) 
iritw) grows algebraically as 

£t(£JcxC (t) for t^<r eq , (20) 

with a small, temperature-dependent exponent a{T) (note that lim^^oo £t(^) = ^(tw), 
so that again). 

The remanent magnetiziation M rem (t) = C(t,0) decreases with increasing domain-size. 
If one assumes M rem (t) oc £t(^) -A ) as for instance in p5|] , one obtains with (|20"D an algebraic 
decay M rem (t) oc t~ a ^ as long as £t(^) "C £,rmeq(T), which would concur with the prediction 
li~6| . However, we find that the time dependence of the remanent magnetization can be better 



fitted to a stretched exponential, similar to the lower bound given in in [T7 



Starting from a one-dimensional domain growth model the exact scaling function for the 
two-point correlations has been predicted in [|1^] . It is (for an unbounded bond-distribution) 
G(r,t w ) = c(r / LT(t w )) (similar to (|9|)) with 

c{x) = [ dy 2y 2 exp(-2x/y) log(y/(l - y)) (21) 
Jo 
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and L/T(t w ) the average domain-size after a waiting time t w at temperature T. Note that 
c(x) ~ exp(— 2.19x) for x G [0,10] to great precision. By rescaling the variable x in fl2~T|) 
by a factor of roughly one half one gets a good agreement with our data, as is shown in 
figure |5|. This means that the free parameter Lq in (Tj| , which is identified with the average 
domain-size, is approximately twice as large as our correlation-length Lq « 2£t(^)- Note 
that due to the construction of the model [113 no prediction about the time dependence of 
this length scale could be made. 



IV. DOMAIN-STATISTICS 

The slow domain growth in spin glasses in general p5| , ^7| and in the random bond Ising 
chain in particular |16jT?|,20| is the reason for their glassy dynamics. This fact has already 



been established numerically for instance in the site-diluted Ising model ||28j| , the random 
field Ising model [^9| and the random bond ferromagnetic Ising model |3(| in two dimensions. 
In contrast to the ±J chain, where this slow dynamics is appropriately described by simple 
kink-diffusion and -annihilation |22[] (which is also the ultimate reason for the " random- walk- 
result" (]T0|) ) , the quenched randomness gives rise to a huge number of metastable states and 
(free) energy barriers of all sizes. Thus the mechanism of domain growth in the presence 
of disorder is significantly modified and slowed down. However, it is still describable by 
means of diffusion and annihilation of kinks (broken bonds) in a random potential |16| . 



Nevertheless it has not been possible up to now to derive an analytically exact expression 
for the average domain-size as a function of the waiting time. 

We calculated the time- dependent domain-size distribution during our Monte-Carlo sim- 
ulations in the following way (note that this definition of domains is different from the 
cluster-definition given in |31j]) : By recording the time averaged local magnetizations 

-r 2t w — 1 

mi(t w ) = - E Si(t) ( 22 ) 
tw t=t w 

for each sample we identified domains as connected segments of the chain in which all local 
magnetizations have the same sign. The length of one domain is just the size of one segment. 
Then we determine the number rii of domains of length I in all samples (a = J2i n i — N ■ S, 



S 



N =system size, S = number of samples) and get the probability P(l,t w ) for a spin to be 
within a domain of length I via 



P T (l,t w ) = I rii(t w ) I S . 



(23) 



In figure |B| we depict the result for two different temperatures. According to the scaling 



behavior found in the last section the probability distribution for the domain-size should 
scale in a similar way 



An appropriate choice for Lx{t w ) is the mean value or average domain-size L(t w ) = 
J dll Pr(l,t w ), a scaling plot is shown in figure [7| and the insert shows LT{t w ) as a func- 
tion of the waiting time. The shape of the curve is similar to that of £t(£u>) 5 although the 
average domain-size is much larger than the correlation length. This is due to the definition 
of domains that we have chosen: It is only the sign and not the magnitude of the local 
magnetizations that determines whether a spin belongs to a domain or not. If a spin has a 
rather small magnetization during a time interval of length t w (which means that it has been 
dynamically "active" rather than frozen) it contributes only marginally to the correlation 
function, but fully to our domain-size. Therefore both quantities are expected to be related 
to each other, but are not equal. 

It is worth stressing that PT{l,t w ) is not Gaussian and that the average LT{t w ) is dom- 
inated by the long tail of the domain-size distribution. The typical domain-size, which is 
defined by the length I, where Pt(/, t w ) attains its maximum, is much smaller than the 
average Lx(t w ), however, it depends on the waiting time in the same way as the average 
(different from what has been observed in the random field Ising chain in a transverse field 



The number of kinks Nx{t w ) in the system, i.e. bonds between spins that are anti- 
parallel, is connected to the average domain length via L T {t w ) = N/N K (t w ). By inspecting 
the probability distribution Kt{J, t w ) for the kink-strength J one obtains information about 
the thermal energy of the chain. The ground state energy per spin of the system is given 
by Eq = — Jo°° dJ P(J), for non-zero temperatures (and for all waiting times t w ) broken 




(24) 
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bonds will increase the thermal energy by an amount AET{t w ) = / °° dJ J KT(J,t w ). In 
figure H we show the kink-strength distribution f^(J, t w ) for two temperatures. We note 
that, comparing it for T = 0.4 with figure [| the kink-strength distribution seems to reach 
stationarity (i.e. independence of t w ) much faster than the domain-size distribution. Since 
the latter yields information about the magnetization of the system, this is in full agreement 
with 1TB] , where it has already been noted, that the relaxation of the thermal energy is much 
faster than that of the magnetization. The whole kink-strength distribution scales nicely 
with the characteristic energy scale 5E T (t w ): 

A ™-^K^b)' (25) 

as is shown in figure || The insert shows the waiting time dependence of the excess energy 
AEx{t w ) in a log- log plot, indicating an algebraic decay over the range of time-scales con- 
sidered (note that it has to saturate at a non-vanishing value for non-zero temperatures, 
which means that an algebraic decay cannot hold forever). This is in agreement with what 



has been proposed in [16 



V. DISCUSSION 



Before summarizing our results we want to discuss the issue of an overlap length f2Ti| , |2l)| , |2"7|| 
in the random bond Ising chain. It has been found that the spin-glass state is extremely 



sensitive to temperature- or field changes by an infinitesimal amount []33|]. This sensitivity 
is observable via the overlap correlation function 

KT,H,ATAh( r ) — [(SiSi+r)^ h (SiSi +r )'^ ±ATiH±Ah ] aiV 

(26) 

cx exp{-r/X Tih (AT, Ah)} , 

where the length scale A (AT, Ah) is the so called overlap length. For a system with a finite 
correlation-length at temperature T and field-strength h, e.g. two-dimensional spin glasses 
or the random bond Ising chain, this is a quite trivial observation. For instance in the ± J 
chain at h = Ah = (see section 2) it easy to see that 

At(AT) = {[ log tanh [J/T]| + | tanh [J/(T + AT)]]}" 1 . (27) 
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The nontrivial implication of (|2q) is that it is expected to hold also within the spin glass 
phase of a true (e.g. higher-dimensional) spin glass, where the correlation-length is infinite 
fl2li| , |2"T||£l]j . The reason is a full rearrangement of the minima of the free energy on length 
scales larger than A by changing the external parameters like T or h only slightly. In |25| 
such a scenario has been introduced artificially (however, motivated by renormalization 
arguments) into the ±J chain via an explicit, stochastic temperature dependence of the 
bonds Jj. 

In the one-dimensional spin glass model such a rearrangement of the minimum-energy 
spin-configurations can be observed at zero temperature (where the correlation length is 
indeed infinite) for a situation, in which one switches on a small external field of strength 



Ah ||16|| . In the presence of a field the ground state of the chain is fragmented into segments, 
in which all spins are aligned with one of the two possible configurations with minimal 
energy at zero-field. In this process weak bonds have to be broken and field energy can be 



gained and thus a new length scale A emerges. It can be shown [|T6| that the length of these 
segments, or the domain-size, depends on Ah via 

X(Ah) oc Ah 2 / 3 , (28) 

which can be identified with the overlap length fl2"S|) for this particular situation T — 0, h — 0. 



At finite temperatures the situation remains similar for small field changes ||16|| . Furthermore 
a generalization of the static overlap-correlation function to a non-equilibrium situation is 



straightforward (see f20|), but has not been investigated yet - - neither numerically nor 
analytically. 

In conclusion we have presented a detailed numerical study of the non-equilibrium dy- 
namics of the random bond Ising chain. We have shown that many of the aging phenomena 
reported in the literature — as for instance strong memory effects, t/t w scaling of the au- 
tocorrelation function and slow domain growth — can also be found in this non-frustrated 
model and also, although less pronounced, in pure Ising models (see section 2, and also Po| . 
where a similar scaling of the autocorelation function as in ( |TT| ) has been derived within an 
£-expansion around four dimensions in a field theory for the non-equilibrium dynamics of 
ferromagnets). 
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We established simple scaling of the two-point correlation function and showed that the 
corresponding (temperature-independent) scaling function concurs with that proposed very 
recently in a one- dimensional domain-growth model [I8|. In addition we have obtained the 
explicit time dependence of the characteristic length scale. Also the domain-size distribution 
and the kink-strength distribution exhibits simple scaling with an average domain size and 
excess energy, respectively. The time dependence of the latter quantities is algebraic as long 
as equilibration is not achieved, which is very reminiscent of what has been observed in 
higher dimensional spin glass models [ff|,|15|. 

Probably most people will agree that there is an essential difference between the simple, 
albeit disordered, model considered here and other spin glass models — as for instance 
the two- and three-dimensional Edwards-Anderson model and the Sherrington-Kirkpatrick 
model. Apart from the presence of frustration the most significant feature of the latter, as 
compared to the random bond Ising chain, is the occurrence of replica symmetry breaking. It 
seems that this equilibrium feature of the SK-model has also interesting consequences for the 
non-equilibrium dynamics ||[10|. Nothing of this will ever be observed in the random bond 
Ising chain. However, many experimental observations on real spin glasses are reminiscent 
of what we have reported here. Thus, as has already been noted in [OT, it seems to us to be 
very useful to look closer into one- dimensional models to catch some of the physics in three 
dimensions. 
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FIGURES 

FIG. 1. The autocorrelation function C(t,t w ) for the ±J Ising chain at zero temperature, 
obtained from the exact result (5). 

FIG. 2. Autocorrelation function C(t,t w ) at different temperatures. 

FIG. 3. Scaling plot of C(t,t w ) versus t/TT(t w ) at temperature T = 0.1. the insert shows the 
characteristic time-scale Tx{t w ) in a log-log plot 

FIG. 4. Two-point correlation function Gt(t, t w ) for different temperatures (T = 0.2 left and 
T = 0.4 right) 

FIG. 5. Scaling-plot of G T (r,t w ) versus v/^T{tw) f° r T — 0-4- The full line is the theoretical 
prediction scaled with a factor 0.46 (see text). The corresponding scaling plot for T = 0.4 yields 
an identical scaling function. The insert shows the correlation length S,r(t w ) in a log- log plot. 

FIG. 6. Domain-size distribution function Pt{1, t w ) for different temperatures(T = 0.2 left and 
T = 0.4 right). 

FIG. 7. Scaling-plot of Pr{l,t w ) versus L,T(t w ) at T = 0.2. The insert shows the average 
domain-size L,T(t w ) in a log-log plot. The corresponding scaling plot for T = 0.4 yields an identical 
scaling function. 

FIG. 8. Kink-strength distribution KT(J,t w ) for different temperatures (T = 0.2 left and 
T = 0.4 right). 

FIG. 9. Scaling-plot of K T (j,t w ) versus AE T (t w ) for T = 0.2. The insert shows AE T (t w ) in 
dependence of t w in a log-log plot. The corresponding scaling plot for T = 0.4 yields an identical 
scaling function. 
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